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ABSTRACT 


In this thesis, high-energy electron generation by stimulated Raman scattering of 
intense laser light is investigated in computer simulations. These high-energy electrons 
can be used to produce a high-energy x-ray source for Nuclear Weapons Effects Testing. 
The simulation results are compared with actual experiments conducted at the Atomic 
Weapons Establishment in the United Kingdom. The simulations show that the heated 
electron temperature is significantly lowered and the frequency spectrum of the scattered 
light is broadened by strong ion density fluctuations in the plasma. These density 
fluctuations are produced when the Raman scattered light undergoes Brillouin 
backscattering. Interactive Data Language is also used to analyze the spectrum of Raman 
scattered light in some recent experiments with the OMEGA laser at the University of 


Rochester. 
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I. INTRODUCTION 


In the Inertial Confinement Fusion (ICF) program, there has been a remarkable 
development of high power lasers, culminating in the National Ignition Facility now 
under construction at Lawrence Livermore National Laboratories. Such lasers are able to 
heat matter to very high temperatures and enable a great deal of high energy density 
physics to be explored in the laboratory. In this introduction, we will briefly discuss 
inertial confinement fusion and give an overview of this thesis work that considers 
another possible application of the National Ignition Facility (NIF). In the next two 
chapters, plasmas, plasma interactions, waves and instabilities, and plasma simulation 
will be discussed. In subsequent chapters, computer simulations of high-energy electron 
generation by intense laser light in underdense plasmas will be presented and discussed. 
Comparisons will then be made between our computer simulations and some real 
experimental data collected at the Atomic Weapons Establishment (AWE) and _ the 
Laboratory for Laser Energetics (LLE). Analysis of some measured scattered light 


spectra from the OMEGA laser at LLE will also be given. 
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Il. INERTIAL CONFINEMENT FUSION & NWET 


Since the late 1940’s the U.S. and other countries throughout the world have been 
searching for a means to replace our over reliance upon fossil fuels. Some of the 
methods currently being explored include solar voltaic cells, high efficiency wind power 
plants and wave energy converters that, in general, produce on the order of tens of 
kilowatts. An alternative path towards fossil fuel independence is inertial confinement 
fusion (ICF). This method, which can lead to power plants generating 1,000 megawatts 
of energy will be investigated using the National Ignition Facility located at Lawrence 
Livermore National Laboratories. Another important use for the NIF facility will be to 
conduct experiments to support the Stockpile Stewardship and Nuclear Weapons Effects 
Testing (NWET) programs. This thesis work is motivated by the potential use of NIF to 
generate a high-energy x-ray source for NWET. 

There are currently two methods that are under consideration for inertial 
confinement fusion. The first approach, indirect drive, is characterized by directing the 
incoming laser beams onto a small cylinder constructed of some high Z material, such as 
gold. The surface of the cylinder starts to ablate; creating plasma that emits high-energy 
x-rays. The x-rays in turn drive the heating and compression of a deuterium-tritium 
target to fusion conditions. The second approach, direct drive, involves heating and 
compressing the fuel capsule by directly shining the laser beams onto the capsule. Each 
approach has advantages and disadvantages. Two major advantages of a direct drive 
method is the possibility of higher gain than an indirect-drive ICF and the ease at which a 
direct drive fusion power plant could be constructed in comparison with an indirect-drive 
power plant. On the other hand, indirect drive is less sensitive to non-uniformities in the 
laser beams and to the Raleigh-Taylor instability during the compression of the fuel 
capsule (Ref 1). The principal approach to ICF to be investigated with NIF will be 
indirect drive. 

Another possible application for NIF is to produce a high-energy xray source for 
NWET, which motivates this thesis work. The x-rays would be produced by 
bremmstrahlung of high-energy electrons. The high-energy electrons are produced 
through the coupling of intense laser light that propagates through the underdense 


plasma. This nonlinear phenomenon is known as stimulated Raman scattering (SRS). 
a) 


In stimulated Raman scattering the laser light will decay into a scattered light wave plus 
an electron plasma wave. After a period of time, the plasma wave has increased in 
amplitude enough to accelerate electrons to a high velocity, making high-energy 
electrons. This process is most efficient when conducted with the density near 0.25 n,, 
where Ne; is the critical density. In theory, almost half of the laser light can then be 
transferred into high-energy electrons, an effect observed in early experiments using the 
1.06 um Shiva laser irradiating underdense plasmas within a hohlraum target. There is 
significant interest in researching the possibility of obtaining an efficient generation of 
high-energy electrons utilizing the National Ignition Facility. It is designed to operate at 
shorter wavelengths than Shiva (planned to be 0.35 [Um and possibly 0.53 um), but 
simulations have shown great promise in this area for NIF if the light is focused to 
sufficiently high intensity. 

To explore this possibility, the hot electron generation by stimulated Raman 
scattering in plasma rear 0.25 n, is investigated using computer simulations that utilize 
an electromagnetic, relativistic particle code. These computer simulations continue 
earlier thesis work by Major Mike Ortelli, USA and Captain David Jones, USA (Ref 2,3). 
Comparisons are made with experiments using 0.53 [wm light at the Atomic Weapons 
Establishment (AWE) laboratory in the United Kingdom. About 50% absorption into hot 
electrons with an effective temperature of ~25 keV is found, in reasonable agreement 
with experimental results. The simulations show that the heated electron temperature has 
been significantly reduced and the frequency of the scattered light significantly 
broadened by the excitation of strong ion density fluctuations in the non linear state. 
These fluctuations appear to be driven by the stimulated Brillouin scattering of the 
Raman scattered light (Ref 4,5). 

Finally, stimulated Raman scattering is monitored in experiments by measuring 
the wavelength spectrum of the back-scattered light as a function of time. An example of 
such experimental data is analyzed using Interactive Data Language (IDL). The resulting 
spectrum can be used to monitor the density of the plasma and to show how the 


efficiency of the scattering varies with time. 


Il. PLASMA 


Plasma is sometimes called the fourth state of matter. Although it makes up 
over 98% of the matter in the universe, in a strange quirk of physics it does not appear in 
a natural state on Earth. We utilize plasmas everyday, from the energy efficient 
fluorescent lights in your kitchen to the cathode ray tubes in older televisions. You can 
see it nearly daily in the Midwest in the form of lightning strikes, yet in practice it is 
usually only seen inside laboratories. Plasma is a quastneutral gas of charged and 
neutral particles that exhibits a collective behavior. The collective behavior comes from 
the fact that the plasma particles are charged. As these particles move around, they 
generate a local concentration of positive or negative charges, which gives rise to electric 
fields. The motions of the charges also induce magnetic fields. These fields affect the 
motion of charged particles far away. It is the long-ranged Coulomb force, which results 
from the solid angle volume of the affected plasma, which gives the plasma a myriad of 
possible interactions. (Ref 6) 

A requirement that the Coulomb force of a neighboring particle be much less than 
the cumulative long range Coulomb force exerted by many distant particles will be 
satisfied when there are many particles in a Debye sphere. A Debye sphere is a sphere of 
radius one Debye length. A Debye length is a typical length over which a charged 
particle's bare electric field has significant influence. In terms of the density and 


temperature of the plasma components, the Debye length, - can be defined by: 


nes} nd 


where n, qg and T are, respectively, the number density, charge and temperature of the 
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particle species, either electrons or ions; «, is Boltzmann's constant and ¢, is the 
permittivity constant in vacuum. 

For distances exceeding the Debye length the electric field of an individual 
charged particle is effectively shielded out by the surrounding plasma. If ® is the 
electrostatic potential of a charge at rest in a plasma, then it can be shown that: 


q elk 
Aner 


where r is the radial distance from the charge. This potential differs from that of an 
isolated charge in a vacuum only by the exponential term. So, the Debye sphere, which is 
a sphere with radius equal to the Debye length, centered on the charge, is a charged 
particle's typical "sphere of influence” on neighboring plasma particles. Not surprisingly, 
the number of plasma particles that occupy this volume plays an important role in 
defining the Coulomb collisional properties of the plasma and determines the significance 


of discrete particle effects in general (Ref. 6). 


When discussing plasmas it is common to refer to the temperature of the plasma 
in terms of energy due to the close relationship between the two properties. The average 
kinetic energy of a plasma is evaluated by 1/2 kgT per degree of freedom, where k zis the 
Boltzmann’s constant and T is the plasma temperature in Kelvin. It is possible, and quite 
normal, to have several different temperatures within any given plasma. This is caused 
because electrons and ions have different temperatures, T; and T, and thus their 
Maxwellian distributions are different. If there were a magnetic field applied, two 
different temperatures within the same species can occur, since forces act differently 
about the parallel and perpendicular axis. Typical plasma may have temperatures on the 


order of a million Kelvin or approximately 100 eV. 
leV = 11,600 K 


A. PLASMA INTERACTIONS 

The fact that plasma particles behave collectively means that plasmas can support 
a wide variety of wave motions and oscillations. One such basic oscillation arises if a 
group of electrons is slightly displaced from their equilibrium positions. The displaced 
electrons feel an electrostatic force seeking to return them to their equilibrium positions, 
but upon arrival there they now have a kinetic energy equal to the potential energy of 
their initial displacement. The electrons overshoot, reconvert their kinetic energy to 
potential energy and a simple oscillation is set up. The frequency of this fundamental 


oscillation is known as the plasma frequency, O ne and is defined by: 





where n, is the mean plasma density and jm, is the electron mass. (Ref. 6) The 
oscillation is of such a time scale that the ions do not have an opportunity to react to the 
oscillating field and they are assumed to be fixed. The equation above is derived with the 
following assumptions: 1) there is no magnetic field, 2) there are no thermal motions 
(KgT=0), 3) the ions are fixed in space in a uniform distribution, 4) the plasma is infinite 
in extent and 5) the electron motions occur only in one direction. This plasma oscillation 
can most properly be described as a plasma wave, because the fringing electric field 
causes a coupling of the disturbance to adjacent layers, and the oscillation does not stay 
localized. Additionally, when electrons stream into adjacent layers within a plasma, their 
thermal velocities will bring information about what is happening in the oscillating region 


of the plasma. 


The angular frequency and the wave vector are functionally related to one 
another, at least in linear theories of plasma waves, by a dispersion relationship. 
Knowledge of the dispersion characteristics of the propagating waves is certainly 
necessary for an understanding of the plasma state. In a simple one-dimensional problem 


the dispersion relationship is as follows: 
2g? 43 e0,2 
@ =Ore 2 th 


where @ is the wave frequency, k is the wave number, @,,, is the plasma frequency and 


the thermal velocity is defined by, Var =2k,T,/m,. This additional thermal correction is 
known as the Bohm-Gross correction. The frequency is now dependent on k and the 
group velocity is finite. There is an imaginary component from the derivation of the 


equation above which results in a linear damping term known as Landau damping. 


Landau damping is a physical effect that is related to the details of the underlying 
velocity distribution functions of the plasma particles. It occurs in the absence of short- 
range binary collisions and hence it is sometimes called collision-less damping. In 
plasma, an initial disturbance is damped (in the absence of binary Coulomb collisions) as 
it propagates away from its point of origin. This can be explained by the following 
argument. Those particles that interact most strongly with a plasma wave are those whose 
velocities are closest to the wave phase velocity. These particles see an approximately 

f| 


stationary electric field. In a stable thermal distribution, for any particular wave speed, 
there are always more particles traveling slightly slower than the wave than those 
traveling slightly faster than the wave. If the particles traveling slower than the wave are 
accelerated by the wave and take energy from it, and those that travel faster than the 
wave give energy to the wave, then there are more particles taking energy from the wave 
than those giving energy to it and the wave damps. Electrons that approach the phase 


velocity of the plasma wave will take energy out the wave and thus gain momentum (Ref 
4). 


With our basic background in some of the definitive properties of plasma, we can 
now investigate how electromagnetic waves in the form of a laser light interact with 
plasma and create the instability cycle which perpetuates the increase in the magnitude of 


the plasma wave through Raman scattering. 
B. PLASMA INSTABILITIES 


First of all though, we should look at frequency and wave number matching 
conditions required for Raman scattering to occur. Raman scattering can be seen as 
simply the resonant decay of an incident photon into a scattered photon plus an electron 
plasma wave (or plasmon). The frequency and wave number matching conditions then 


become: 
0, = 0, + Oo 


k,=k,+k 


plasma 


where O, and k, are the frequency and wave number for the incident light wave, w, and 


k,are for the scattered light wave, while @,, and k 


plasma te for the electron plasma wave. 
It is clear that the minimum frequency of light within the plasma must be at least Wpe, 
therefore 2 2@pe. Light will not propagate if Wpe > @. Because n/Mer = (Wpe/@o), then 
Ne S 0.25ner , where n, is the plasma density and n, is the critical density (Ref. 4). The 
fraction of energy that is deposited in the plasma wave heats the plasma as the wave 
damps. This kinetic energy increase is what produces the high-energy electrons that are 
characteristic of the Raman instability. However, since these electrons can preheat the 


fuel in laser fusion applications, this instability has to be monitored carefully. 
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The Raman instability provides a mechanism by which the electron plasma wave 
and the associated scattered wave can build up in amplitude utilizing a feedback system. 


Figure | shows the stages of the instability cycle. 


During the first stage any perturbation, such as thermal noise can cause a small 
charge density fluctuation to occur. This density fluctuation, 6n, is due to a thermal level 
plasma oscillation. The oscillatory motion of electrons in the light couples with 6n to 
induce a transverse current fluctuation: 


_ xe E,8n 


m0), 


oJ 


where e is the electron charge, FE j is the laser electric field, m,is the mass of an electron 


and q@_, is the light wave frequency. 


The perturbed current produced, 6J, will then drive the scattered light wave at the 
beat frequency. The scattered light wave grows from noise as a perturbed wave with 
electric field amplitude, 5E (a light wave at an up shifted or downshifted frequency). 
Finally this electric field will beat in turn with the electric field produced by the laser to 


create a perturbation in the electric field pressure. 
E > € 5 
OP field = ae: +6E) ae =¢,E,6E 


This pressure perturbation will then reinforce the initial density fluctuation. This 
feedback system leads to exponential growth of both the plasma wave 6n) and the 


scattered light wave (SE) (Refs. 4,7,8,9). 


én 


8P,_jeia * ESE SJ < E,&n 


SE 


Figure 1. Plasma instability feedback loop for Raman backscatter growing from noise. 


A second important instability is stimulated Brillouin scattering, which involves 
the production of a scattered light wave plus a low frequency ion acoustic wave. The 
production of the ion acoustic wave is similar in manner to the production of the Raman 
instability, except that resonant decay of the incident photon results in a scattered photon 


plus an ion acoustic phonon. The frequency and wave number matching then becomes: 
@, =0,+0,, 
k,=k,+kiq 


where @,, and k,, are for the ion acoustic wave. The frequency of the ion acoustic wave 
is much less than either the plasma frequency or the incident photon frequency. Therefore 
this instability can exist at all densities up to the critical density. Brillouin scattering is a 
significant concern in laser fusion research, because nearly all of the energy can be 
transferred to the scattered light wave. This may appreciably degrade the rate of 


absorption, and SBS can also lead to deep density modulations in the plasmas that will be 
10 


discussed later. If the wavelength of the ion acoustic waves is considered to be small by 


taking k°A;, >>1, then the dispersion relationship becomes: 


2 
2 noe 


oO 


eM 





where M is the ion mass. There becomes a complementary behavior between the ion 
acoustic waves and the electron plasma wave. At higher wave numbers the ion acoustic 
waves exhibit a constant frequency, but become constant velocity at smaller k values. 
The electron plasma waves are constant velocity at high k values, but become constant 
frequency at small k values (Ref 4,5,6). This is shown graphically in Figure 2, where v,, 
is the thermal velocity, v,is the ion acoustic velocity and , represents the ion plasma 


frequency. 


Electron 





Figure 2. Comparison of the dispersion curves for electron plasma and ion acoustic waves 
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IV. GREEN LIGHT COMPUTER SIMULATIONS 


The advancement in the past three decades in laser power output and computer- 
processing speed has played a critical part in the ability to analyze laser-plasma 
interactions. As an example, the development of the SHIVA laser by Lawrence 
Livermore National Laboratories (LLNL) in 1978 was one of the first steps in creating a 
mega-joule class laser. SHIVA was rated at 10 kilojoules and was superceded by the 
NOVA laser, which was rated at 40+ kilojoules, in 1984. Today construction continues 
on the NIF laser at LLNL, which will be used for laser fusion research and stockpile 
stewardship. NIF is expected to be fully operational by 2008 and will have a 1.8 mega- 
joule output of 0.35 um light. This means that on a rough average, the laser capacity 
output of some of the largest lasers in the world doubles about every three years. Similar 
in comparison is the progression of computer processing. Gordon Moore, co-founder of 
Intel made an observation in 1965 that the number of transistors per square inch on 
integrated circuits had doubled every year since the integrated circuit was invented. 
Moore predicted that this trend would continue for the foreseeable future. In subsequent 
years, the pace slowed down, but data density has doubled approximately every 18 
months. This is the current definition of Moore's Law. Most experts expect Moore's Law 


to hold for at least another two decades. 


Computer simulations have become an extremely important facet in the analysis 
of laser plasma interactions. In the early years of computer simulation of plasmas, the 
“supercomputers” at the time could only handle on the order of a thousand or so particles. 
These room-sized computers could only manipulate the most rudimentary calculations 
and took a great deal of time to do so. Today, an ordinary laptop or desktop computer 
with a sizeable amount of memory can easily handle hundreds of thousands of particles 
and larger mainframes can handle very complex multidimensional problems with 


millions of particles. 
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The Department of Energy (DOE) is currently utilizing ASCI White, a IBM 
constructed computer at LLNL capable of calculating 12.3 trillion floating point 
operations per second (Teraflops), to simulate nuclear weapons testing. These 
simulations are essential due to the U.S. observation of the Comprehensive Test Ban 
Treaty. Without the possibility of conducting underground tests, the computer 
simulations help make sure that the nuclear weapons stockpile is safe, reliable and 
operational. With ASCI White, we can analyze the factors involved in a nuclear 
detonation. In 2004, there will be a 60 Teraflops supercomputer known as ASCI Purple 
at LLNL. The current ASCI project is expected to have a 100 Teraflops supercomputer 
by about 2006. As a comparison of the progress in computing power, in 1985 it would 
have taken the fastest supercomputer 60,000 years to perform the full three-dimensional 
nuclear simulation that this future supercomputer will be able to accomplish in less than a 


month! 
A. Particle in Cell Code Simulations 


Computer simulations using Particle in Cell (PIC) codes offers scientists a 
powerful and direct approach in the investigation of kinetic and non-linear effects. An 


example of the basic cycle of the PIC code is shown in Figure 3. 


{7} 






Newton’s Laws Grid Assignment 


{EB} 7} 


Maxwell’s Equations 


Figure 3. _ Particle code simulation cycle. 


It would be constructive to cover the components of a very basic PIC code. For 
simplicity, we treat ions as a fixed and neutral background in this introductory discussion. 
We can sometimes make this assumption because the mass of the ions is so much larger 
than that of the electrons. Secondly, we impose no magnetic fields. In the limit of 
electrostatics, the magnetic field generated by plasma currents is negligible. Hence, we 
only consider electrostatic fields and Maxwell’s equations reduce to e,Ve E =o, where 6 
is the charge density. In the interest of producing a very simple code, variations are 
allowed only in one direction. The problem starts off by calculating the charge density. 
This is accomplished by dividing the system created into a number of cells that are 
determined by the user in the computer code. The cell size, 5, is then made unity which 
allows L=N., where L is the length of the system and N, is the number of cells within 
the system. Periodic boundary conditions will be used in this simplest example, meaning 
that charges leaving the system on the left (right) boundary are simply reintroduced on the 
right (left) boundary. The charge can simply be assigned to its nearest grid point location; 
however it is better to share the charge between its two nearest grid points (Figure 4). For 


a charge located a distance Ax to the right of the {” grid point, you have: 
Ao(i) = gd —-Ax) 
Ao(i+1) = qAx 
where q is the basic electron charge and Ao is the increment in the charge density. Once 


the charge is assigned on the spatial grid, the electric field is determined next using 


Poisson’s equation in the one-dimensional form. 


£, = =n, =N,) 


where n, is the ion density and n, is the electron density. 
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One possible method is to Fourier transform the charge density and Poisson’s 
equation. The inversion of the Fourier transform will give the electric field on the spatial 


grid. An easier method is to use finite differences. 


o(i+)+o(i 
EG+D=E@OL Fe pace 
The electrical force is assigned to each of the particles using the same method 
used to assign the charge to the grid. By considering a particle a distance Ax to right of 


the i grid point, the force F on the particle is determined by linear interpolation. 
F=qE(@[1—-Ax]+ gE(i+ 1lAx 


The velocity, v, and position, x, of each particle are then advanced At in time using a 


“leap frog” algorithm. 


aD = a + F"At 
1 
x™ ax" 4p" 2A 
Each of the superscripts denotes a time step. If x and v are defined as being one-half time 
step apart, then second order accuracy is achieved. Initially, position and velocity can be 
defined at the same time, but the velocity can be displaced backwards using the force at 
t=0. 

The plasma cycle continues to be computed by using a time step small enough to 
resolve the characteristic oscillations of the plasma, usually the electron plasma 
frequency. The time step commonly used is on the order of 0.2@,-. A grid size of about 
Apex helps resolve the scale lengths characteristic of the collective behavior. Aliasing is a 
numerical instability that arises when the wave number k cannot be distinguished from 
spurious ones with wave numbers k+(27n/5), where 6 is the grid spacing and n is an 


integer. This can be eliminated by not choosing the grid size to be a small fraction of the 


Debye length (Ref 10). 
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Figure 4. Particle in cell code sharing method 


B. TURBOWAVE AND THE AWE EXPERIMENTAL ANALYSIS 


One of the possible uses for the National Ignition Facility (NIF) is to generate 
high-energy electrons with an effective temperature of between 20-50 keV using either 
0.53pm or 0.35pm laser light in order to produce a high-energy x-ray source for nuclear 
weapons effects testing. Computer simulations are used to investigate generation of these 
high-energy electrons via stimulated Raman scattering (SRS) of intense laser light (Ref 
8,9). It is shown that a fairly efficient hot electron generation (~50%) can be achieved if 
plasma with a density near 0.25 n, is irradiated with sufficiently intense laser light. The 
calculations compare favorably with the hot electron generation observed in recent 
Atomic Weapons Establishment (AWE) experiments in the United Kingdom using 


intense 0.53 um light. Similar results are theoretically expected for 0.35 tm light 
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provided that the plasma density and laser intensity are scaled about a factor of two 


higher. 


The computer simulations are carried out using Turbowave, a parallel 3-D, 
electromagnetic, relativistic particle code written by Dan Gordon from the Naval 
Research Laboratory (NRL). In contrast to the simple case discussed to illustrate the 
concepts, this particle in cell (PIC) code solves Maxwell’s equations along with the 
relativistic equations of motion of many electrons and ions, which represent the plasma. 
One and two-dimensional simulations were used to follow the propagation of intense 
laser light through a slab of plasma, as indicated in Figure 5. Periodic boundary 
conditions are used in the other directions. Particles striking the left or right boundaries 
are re-emitted with their original thermal distribution. The laser light is sent in from the 
left boundary, and light waves striking the left or right boundaries are allowed to pass out 


of the system (Ref 11). 
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Figure 5. _ Intense laser light interactions with plasma slab 
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1. Atomic Weapons Establishment HELEN Experiments 


Many simulations have been carried out with a wide variance in plasma densities, 
laser intensities, and electron and ion temperatures. We will discuss an instructive 1-D 
simulation that was motivated by the recent experiments at AWE. For this simulation, 
the initial uniform plasma density is 0.23 ng, where ng is the density at which the 
incident light striking the plasma is totally reflected. We are choosing a uniform, 


collisionless plasma with density @,, /o, =1/2.1=.Jn,/n, . The following other 


parameters are also considered: 


e The electron temperature, Tz, is 2.5 keV, the ion temperature, Tj is 0.625 
keV, 
e The size of the uniform plasma slab is 68.5 free space wavelengths for the 


laser light or a plasma length of 204.8c/o,,, 


e v,,/¢ =0.04 (the driver’s amplitude, related to the dimensionless laser 


intensity; Uos is the electron’s velocity as a result of the laser’s electric 


field) 


e Duration of the run is 3277/@,, OF 1.94 x 10°” seconds, which is 


sufficient for the rapidly growing instability to saturate. 


As the laser light propagates through the plasma slab, it undergoes stimulated 
Raman scattering, i.e., it decays into a backscattered light wave plus an electron plasma 
wave with a frequency that is principally @pe, the electron plasma frequency, plus a small 
thermal correction. Due to the instability discussed earlier, the electron plasma wave and 
the scattered light wave will continue to grow in amplitude. The fraction of energy that is 
deposited into the wave is responsible for the increase in kinetic energy that results in the 


production of hot electrons. 


The efficiency of the electron heating can be obtained by monitoring the various 


energy histories in the simulation. Figure 6 shows the evolution of the kinetic energy of 
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the plasma, the energy carried off when the heated particles exit the boundaries, the 
energy in the electric and magnetic felds in the plasma, as well as the Poynting flux 
through the left and right boundaries. At the end of the run, the kinetic energy is about 
0.4 (dimensionless units used), which is up from about 0.15 initially. The lost particle 
energy equates to about 0.3. If these energies are added together and divided by the total 
energy input (the amount of energy crossing the right boundary as if there were no 
plasma present) from the laser, you obtain an estimated absorption. The energy output 
from the right boundary is about 1.15; therefore, the fraction of the incident laser energy 


absorbed is about 48.7%. 


—e— KineticEnergy 
—s -LostParides 
— @ -AReldEnergy 
| --»- - LeftPoyntngAux 
- -+- - RightPoyntngAux 


Relative Energies 





ut) 500 1000 1500 2000 2500 3000 3500 


Figure 6. Energy diagram for 2.5 keV computer simulations. 
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For this simulation, about 50% of the incident laser energy is absorbed into heated 
electrons. This compares reasonably well with the analogous experiment, where about 
20% absorption into high-energy electrons was measured. These experiments with the 
HELEN laser in the United Kingdom at AWE were conducted utilizing various average 
plasma densities and laser intensities irradiated with a single 20 unsmoothed beam (Ref. 
10). Figure 7 shows the fraction of the laser energy absorbed into high-energy electrons 
(faor) as a function of average plasma density. It can be seen that fyor is a function of 
both plasma density and laser intensity, with a peak being centered at about 0.25 no 
where the maximum amount of absorption takes place. Of course, there is actually a 
range of densities in the experiment, so some absorption into hot electrons occurs even 


when the average density exceeds 0.25 n,, the maximum density for which SRS occurs. 
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Figure 7. High fyor only with high intensity, unsmoothed beams 
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The lesser absorption is to be expected as the laser beam had to first form the 
plasma in this particular experiment. We estimate that the ideal conditions for SRS were 
only obtained during the second half of the laser pulse. In addition, the intensity of the 
laser light was somewhat less in the experiment. The equation for intensity for this 


simulation is: 


l= 1.35x10'%a? _ 1.35x10'*(0.04)? 


= 7.63x10'°W/em? 
ic (0.532) 


where a=v,,/c. 


When the electric field of the plasma wave becomes sufficiently large, electrons 
are accelerated to high energies, producing a high-energy tail on the electron energy 
distribution. An example of this high-energy tail is shown in Figures 8a and 8b. In the 
plots you can see the probability distribution of electrons plotted, as a function of energy 
expressed as a percentage of the rest mass energy. The slope of the distribution gives the 
heated electron temperature. For comparison, we look at the high-energy tail without ion 
motion not allowed (essentially making the ions of infinite mass). In Figure 8a, the effect 
of ion motion not being allowed is quite apparent. The hot electron temperature, 


determined from the slope of the distribution at high energy is almost 130 keV! 


In a more realistic scenario, however, we must include ion motion and allow ion 
density fluctuations to develop. Figure 8b shows the results when ions are allowed to 
move. In this case, the hot electron temperature is about 26 keV. This means an overall 
reduction in the hot electron temperature of almost 80%, showing that the ion density 
fluctuations created severely reduce the characteristic temperature of the high-energy 


electrons. 


The hot electron temperature compares very favorably with actual data obtained 
from testing at AWE. The measured Tho during HELEN experiments was approximately 


30 keV. It is important to note that this heated electron temperature has been 
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significantly reduced by the development of pronounced ion density fluctuations in the 


plasma. 


Figure 8a. High energy tail distribution with no ion motion (~130 keV) 
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Figure 8b. High energy tail distribution with ion motion allowed (~26 keV) 
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Figure 8. Hot electron energy distribution for 2.5 keV computer simulation 


The evolution of the plasma density fluctuations is shown in Figures 9 and 10. 
Figure 9a shows the simulation early on when the density fluctuations have not yet 
developed. In Figure 9b, you can see the ion density fluctuations beginning to propagate 
within the first half of the plasma. Figure 10a shows the initial development of deeply 
modulated ion fluctuations. Finally 10b, on a slightly different scale to emphasize the 


deep canyons in the density, shows the final evolution of the plasma modulations. 


The fluctuations reduce the heated temperature since they couple the Raman- 
generated plasma waves into plasma waves with a slower phase velocity, which in turn 
accelerate electrons to lower energies. These ion density fluctuations have a rather long 
wavelength (about three times A», the incident laser wavelength). The fluctuations have a 
simply physical explanation. In the near 0.25 n, plasma, the Raman scattered light wave 
has a frequency near Wpe and a long wavelength (about 6A,). This scattered light wave 
becomes large in amplitude and undergoes simulated Brillouin scattering (SBS), which 
corresponds to its decay into a backscattered light wave plus an ion acoustic wave. The 
ion density fluctuations associated with this ion acoustic wave have a wavelength of 


about 3A, in agreement with the results in Figure 10b. 
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Figure 9. _—_ Plasma density fluctuations for 2.5 keV plasma run (early) 
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Figure 10. Plasma density fluctuations for 2.5 keV plasma run (late) 
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Let’s now consider a simulation with a lower electron temperature, as this may be 
closer to actual experimental temperatures. For this simulation, the initial uniform 
plasma density is again 0.23 n,, where ng is the density at which all incident light 
striking the plasma is totally reflected. We are choosing a uniform, collisionless plasma 


with density ,, Io, =1/2.1=.]n,/n,, . The following other parameters are also 


considered: 

° The electron temperature, T., is 1 keV, the ion temperature, Tj is 0.25 keV, 

e The size of the uniform plasma slab is 68.5 free space wavelengths for the 
laser light or a plasma length of 204.8c/@,,,, 

° v,,/c =0.04 (the driver’s amplitude, related to the dimensionless laser 
intensity; Vos is the electron’s velocity as a result of the laser’s electric 
field), 

° Duration of the run is 3277/@,, or 1.94 x 10°’? seconds. 


With a reduced background electron temperature we find that the amount of 
Raman scattering and the high-energy electron temperature are about the same as before. 
We use the same data analysis as for the earlier computer simulation. The efficiency of 
the electron heating is obtained by analyzing the energy diagram for the 1 keV 
simulation. Figure 10 shows the evolution of the kinetic energy of the plasma, the energy 
carried off when the heated particles exit the boundaries, the energy in the electric and 
magnetic fields in the plasma, as well as the Poynting flux through the left and right 
boundaries. At the end of the run, the kinetic energy is about 0.24 (dimensionless units 
used), which is up from about 0.04 initially. The lost particle energy equates to about 
0.24. If these energies are added together and divided by the total energy input (the 
amount of energy crossing the right boundary as if there were no plasma present) from 
the laser, you obtain the total absorption. The energy output from the right boundary is 
about 1.0, therefore the amount of laser light absorption is about 44.0%, about the same 


as in the previous example. 
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Figure 11. Energy diagram for 1 keV computer simulation. 
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Next we examine the high-energy tail created by the Raman instability. We 
notice that the average temperature of the hot electrons is similar to the 2.5 keV 
simulation. Figure 12 shows the heated electron distribution function. From the slope, 
the heated electron temperature can again be found. For this 1keV input, you receive an 


average electron temperature of 20.4 keV, a bit lower than before. 


High energy tail distribution (~20 keV) 
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Figure 12. Hot electron energy distribution for 1 keV computer simulation 
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Displays of the plasma fluctuations for the 1 keV simulation are shown in Figures 
13 and 14. Figure 13a shows the timeframe in which the plasma fluctuations have not yet 
developed. In Figure 13b, you see the magnitude of the ion density fluctuations starting 
to increase strongly, on a relative scale. Figure 14a shows the development of the 
modulated fluctuations throughout the plasma, somewhat different from the 2.5 keV case 
in which the fluctuations were only deeply modulated in the first half of the plasma. 
Finally Figure 14b shows that the deep modulations are occurring throughout the entire 
length of the plasma. The fluctuations reduce the heated temperature since they couple 
the Raman generated plasma waves into plasma waves with a slower phase velocity, 


which in turn accelerate electrons to lower energies. 
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Figure 13a 
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Figure 13. Plasma density fluctuations for | keV plasma run (early) 
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Figure 14. Plasma density fluctuations for 1 keV plasma run (late) 
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Figure 15. 2. OMEGA Laser Experiments at 15% critical density 


Particle in cell codes, such as Turbowave, are flexible tools to explore many 
possibilities and parameters. Many scientists are using PIC codes to investigate how to 
limit the creation of hot electrons and Raman scattering. Hot electrons help preheat the 
fuel, which then prevents the full compression of the fuel to occur. This is a significant 
problem if you are trying to make laser fusion successful. Others are seeking ways to 
increase the generation of hot electrons in order to conduct Nuclear Weapons Effects 
Testing (NWET) experiments. More hot electrons mean that more high-energy xrays 
will be created from the bremmstrahlung of the electrons with a high Z material. It 
would help both camps if we can understand when SRS is efficient and when it is not. 
Indeed, some recent experiments with the OMEGA laser explored hot electron generation 
for lower intensities and shorter wavelength laser light (Ref. 13). In these experiments, 
gasbag plasmas with various densities were irradiated with 0.35um laser light at an 
intensity of about 10'° W/cm’. Typically fizor was measured to be only a few percent. 

Some simulations relating to these Omega experiments were carried out using the 
Turbowave program. Effects such as much reduced absorption into high-energy 
electrons was readily seen in simulations with lower intensity laser light. However, our 
simulated plasmas were too small and the hot electron generation too weak to provide 
reliable comparisons with experiments. PIC simulations are ideal for investigating 
strongly nonlinear regimes, such as in the AWE experiments at high intensity. It is much 
harder to accurately compute the hot electron generation when it is weak, especially with 
simulations on a MAC G4. One would need to use a supercomputer that would enable 
simulation of very large, weakly driven plasmas over much longer times with many more 
particles. Details such as density gradients in the experiments would also need to be 
considered. Such simulations are beyond the scope of this thesis, which is more focused 
on high intensity experiments where the hot electron generation is strong, as desired for 


an x-ray source. 
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C. ANALYSIS OF THE FOURIER TRANSFORMS OF SELECTED RUNS 


Another useful tool to analyze the computer simulations is to take a Fourier 
transform of the transverse electric field near the left- most boundary of plasma. This will 
show the frequency components that make up the electric field. We will look briefly at 
the previous example; the .23 n,; run that was conducted with an electron temperature of 


2.5 keV. 


Figure 15 shows the results for the n = 0.23nm, simulation in which we would 
expect two peaks, indicating the incident laser light field and the Ramarscattered light. 
In the figure, spectral density (in relative terms) is plotted against o/o pe: AS expected, 
we see a well-defined peak at approximately w/w pe = 2-1, which shows the incident laser 
light. We also see a range of lower frequencies that shows the Ramarrscattered light. If 
we look closer in Figure 16, we see that the range of frequencies is from about 
o/@ pe = -95—1.09. This spectrum is much broader than expected based on linear theory. 
However, such a range of frequencies seem quite reasonable when account is taken of the 
large density fluctuations present in the final nonlinear state in the plasma. Figure 17 
shows the same simulation except ion motion is not allowed. Note the very definitive 
peaks at 2.1, which indicates the incident light, and at about 1.09, which is the Raman 
backscattered light. Clearly the ion density fluctuations cause the broadening of the 


scattered light frequency. 
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Figure 16. Fourier transform of transverse electric field for 2.5 keV 
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Fourier transform localized showing Ramarscattered frequency 
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Figure 18. 
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D. SPECTRAL ANALYSIS OF OMEGA SHOTS UTILIZING IDL 


Finally, we briefly consider an analysis of some experimental measurements of 
the scattered light. This data gives the spectrum and intensity of the back-scattered light 
as a function of time. Such data gives information regarding the density of the plasma 
from which the scattering is occurring. We use Interactive Data Language (IDL) 
programming to analyze laser shots conducted by Robert Kirkwood from LLNL in July 
2002 on the OMEGA laser at the University of Rochester. Major Mike Ortelli from 
DOE, Professor Chris Olsen from the Naval Postgraduate School and myself developed 
the IDL code used. The IDL code takes hierarchal data format (.hdf) files that are 
produced from the online target diagnostic Full Aperture Backscatter System (FABS) and 


converts them into portable network graphics (.png) as seen in Figure 18. 


oovese00 


Figure 21. Stimulated Raman Scattered streak in .png format 





These graphics clearly show the SRS streak as a function of time, with 
wavelength increasing with ascension up the y-axis. You can see a characteristic streak 
at about 532 nm, which shows the incident green laser light from the conversion crystal. 
You can estimate the time scale by looking at the fiducial at the bottom of the figure. 
The distance between each of the marker is approximately 1/2 nanosecond. Therefore we 
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can estimate that the overall length of this laser shot was approximately 1 — 1.5 
nanoseconds. Once we have this image, we can do a spectral density image using the 
IDL program. From this spectral image density we can see the intensity of the scattered 
light as a function of wavelength and time. This image is not exactly correct, as it does 
not take into account the dispersion within the fiber optics of the system. Dispersion is 
caused when different wavelengths of light travel at different speeds within a medium 
(other than vacuum). In the fiber red light travels faster than blue light and this must be 
taken into account. Within the IDL program, we incorporated a modified fiber dispersion 
correction that we received from Wolf Seka (Ref 14), who has conducted similar analysis 
on other OMEGA experiments. This fiber correction is shown in Figure 19, which shows 
the shift in time when dispersion is taken into account for wavelengths of less than 
550nm. In Figure 20, we can see the spectral density plot after the fiber correction is 
taken into account. This plot shows a contour map of the relative intensities. If you 
compare Figures 18 and 20, you can see that plot follows very closely the actual spectral 
density plot downloaded from OMEGA; however in this case the fiber correction did mt 
make a significant difference. When available, an improved fiber correction can be 
incorporated for future data analysis. Finally, in Figure 21 we plot the intensity 
integrated over time as a function of wavelength. In this data, the scattered light peaks at 


a wavelength of about 650 nm, which corresponds to SRS from a density of about .16 n,,. 
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V. CONCLUSIONS 


Stimulated Raman scattering (SRS) represents the resonant decay of an intense 
light wave into an electron plasma wave plus a scattered light wave. When the electron 
plasma wave becomes large in amplitude, it accelerates electrons to a high energy. These 
high energy electrons can be used to produce a high energy x-ray source for NWET. 
Generation of such a x-ray source is another possible application for the National Ignition 


Facility. 


Computer simulations were used to investigate the high-energy electron 
generation by SRS. Efficient generation was found when plasma with a density of about 
0.25 ng is irradiated with intense 0.53 um laser light. The simulations compare favorably 
with experimental data from the HELEN laser in the United Kingdom. In both 


experiments and simulations the heated electron temperature was about 30 keV. 


The simulations show that the temperature of the heated electrons is strongly 
reduced and the frequency spectrum of the scattered light is broadened by large ion 
density fluctuations that form in the plasma. These long wavelength density fluctuations 
are associated with the ion acoustic wave which forms when the Ramarscattered light 
waves undergo stimulated Brillouin scattering; i.e. decay into a backscattered light wave 
plus an ion acoustic wave. In related work, the spectrum of the Ramarscattered light 


measured in some OMEGA experiments was analyzed using an IDL code. 


Finally, the computer simulations also show that the efficiency of high energy 
electron generation is much reduced when lower intensity, shorter wavelength laser light 
is used. This trend agrees with some experiments with the OMEGA laser, where the 
absorption into hot electrons was observed to be only a few percent. Realistic simulation 
of this weak regime would require use of super computers, which could be a future thesis 


topic. 
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